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After  polytetrafluoroethylene  (PTFE)  treatment,  due  to  the  coexistence  of  hydrophilic  and  hydrophobic 
pores,  the  gas  diffusion  layer  (GDL)  of  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  shows  a  character 
of  mixed  wettability,  which  in  turn  affects  liquid  water  and  mass  transfer.  A  series  of  fractal  models  are 
developed  in  this  work  to  investigate  the  effect  of  GDL’s  wettability  on  liquid  water  and  gas  permeation  in 
GDL.  Compared  to  the  widely-used  empirical  models,  the  proposed  model  of  saturation  versus  capillary 
pressure  in  a  good  agreement  with  experimental  data  is  more  suitable  for  the  GDL  of  mixed  wettability. 
By  using  this  model,  liquid  water  saturation  is  found  to  be  positively  correlated  with  tortuosity  and  pore 
area  fractal  dimensions  for  a  hydrophilic  case  and  hydrophilic  pore  fraction,  whereas  to  be  negatively 
correlated  with  these  fractal  dimensions  for  a  hydrophobic  case.  Furthermore,  theoretical  predictions 
on  gas  and  water  relative  permeability  are  made  via  the  proposed  models.  Water  relative  permeability 
increases  with  the  increases  in  the  fractal  dimensions  for  a  hydrophobic  GDL  and  liquid  water  saturation, 
whereas  it  decreases  with  the  increases  in  the  fractal  dimensions  for  a  hydrophilic  GDL  and  hydrophilic 
pore  fraction.  For  gas  phase,  an  opposite  result  is  obtained. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Gas  diffusion  layer  (GDL)  is  a  component  crucial  for  the  opera¬ 
tion  of  polymer  electrolyte  membrane  fuel  cell  (PEMFC),  as  it  plays 
an  important  role  in  the  delivery  of  reactants  from  flow  channels 
to  catalyst  layers  and  the  effective  removal  of  product  water  from 
the  electrode.  At  higher  current  densities,  a  dramatic  increase  in 
reaction  rates  leads  to  the  corresponding  increase  in  water  gen¬ 
eration.  Liquid  water  may  fill  the  pore  network  within  GDL  and 
decrease  the  effective  diffusivity  of  reactants  through  the  layer.  In 
the  actual  application  of  PEMFC,  to  facilitate  the  removal  of  liq¬ 
uid  water,  the  GDL  is  usually  treated  with  a  non-wetting  polymer 
such  as  polytetrafluoroethylene  (PTFE)  to  create  hydrophobic  sur¬ 
faces  and  pores  throughout  GDL,  resulting  in  coexistence  of  the 
hydrophilic  and  hydrophobic  pores  and  change  in  the  wettability 
of  diffusion  medium. 

The  presence  of  excessive  liquid  water  in  the  porous  GDL  blocks 
part  of  pores  towards  the  catalyst,  debilitating  the  performance 
of  device.  Pore  scale  phenomena  associated  with  the  presence  of 
liquid  water  inside  porous  media  and  the  effects  of  liquid  water  on 
capillary  pressure,  relative  permeability  of  the  wetting  and  non¬ 
wetting  phases  have  become  increasingly  recognized  in  the  recent 
years.  Many  experimental  measurements  and  theoretical  analysis 


*  Corresponding  author.  Tel.:  +86  18672952831;  fax:  +86  027  87859589. 
E-mail  address:  a_laly@163.com  (Y.  Shi). 

0378-7753 /$  -  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 2.02.092 


have  been  taken  to  investigate  the  properties  of  porous  GDL  in  the 
PEMFC 

Some  of  the  experiments  have  been  carried  out  to  study  the 
morphologic  properties  of  the  GDL  by  using  mercury  intrusion 
porosimetry  (MIP)  [1-6],  which  cannot  distinguish  between  the 
hydrophilic  and  hydrophobic  components  of  the  diffusion  medium. 
And  most  of  these  studies  have  focused  on  the  dependence  of  capil¬ 
lary  pressure  on  saturation  for  mercury-air  system,  while  few  have 
been  concerned  with  that  for  water-air-GDL  system,  as  required 
in  mass  transfer  model.  Gostick  et  al.  [7]  measured  this  depend¬ 
ence  by  using  standard  porosimetry  method  with  water-air  as 
working  fluid  pairs.  Their  following  work  [8]  exploited  an  appa¬ 
ratus  to  survey  the  relationship  between  the  air-water  capillary 
pressure  and  water  saturation,  which  enables  the  determination  of 
the  hydrophilic  pore  size  distribution.  Experimental  measurements 
contribute  a  lot  towards  understanding  the  capillary  behavior  of 
water-air-GDL  system,  and  always  go  with  polynomial  fits  or 
empirical  formulae. 

For  the  GDLs,  some  researchers  used  a  function  that  they  devel¬ 
oped  from  fitting  data  to  relate  the  capillary  pressure  to  the 
saturation  [9,10,2],  others  adopted  the  van  Genuchten  model  and 
the  Brooks-Corey  model  [7],  while  majority  described  the  relation¬ 
ship  with  the  widely  used  and  empirically  determined  Leverett 
J-function  and  its  modified  formulae  [11-16],  whose  validity  for 
their  use  in  PEMFC  modeling  has  recently  received  strong  critics 
since  they  ignore  the  detail  pore  morphology  of  the  GDL  In  sum¬ 
mary,  there  is  much  debate  on  the  empirical  expression  between 
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Fig.  1.  Fractal  bundle-of-hydrophilic-and-hydrophobic-capillaries  model  for  porous  GDL. 


the  capillary  pressure  and  the  saturation.  Zamel  and  Li  [17]  and  Wu 
et  al.  [18]  compared  several  expressions  to  find  out  their  effects  on 
capillary  pressure,  and  their  research  also  showed  that  the  choice 
of  the  empirical  expression  for  capillary  pressure  is  crucial  for  the 
simulation  on  liquid  water  transport  in  the  PEMFC. 

In  general,  the  relation  between  capillary  pressure  and  water 
saturation  is  used  in  numerical  simulation  for  the  theoretical 
research  on  the  spatial  distribution  of  liquid  water  and  saturation 
in  the  GDLs.  Wang  and  Nguyen  [19],  Zamel  et  al.  [20]  and  Sinha 
and  Wang  [21  ]  investigated  the  effect  of  capillary  properties  of  the 
GDL  on  the  liquid  water  transport  and  water  saturation  level  using 
the  numerical  models.  Most  of  the  theoretical  analyses  employ 
the  empirical  expressions,  although  a  short-coming  is  that  their 
coefficients  do  not  have  a  direct  correlation  to  physical  quantities, 
hindering  the  fundamental  understanding  of  water  percolation  and 
movement  within  the  diffusion  medium. 

Except  for  the  above-mentioned  drawback,  the  models  for  theo¬ 
retical  analysis  mostly  treat  the  GDL  as  entirely  hydrophilic,  which 
means  that  the  liquid  pressure  must  always  be  below  the  gas 
pressure.  It  is  taken  into  account  in  the  models  by  assuming  the  sat¬ 
uration  with  a  value  of  zero  at  the  interface  of  a  diffusion  medium 
with  a  gas  channel,  corresponding  to  the  liquid  pressure  much 
lower  than  the  gas  pressure  (and  even  approaching  a  value  of  zero) 
at  this  interface  [22].  Since  the  GDLs  have  added  PTFE  to  keep  them 
from  flooding,  the  coexistence  of  hydrophilicity  and  hydrophobic- 
ity  of  GDL  results  in  a  mixed-wetting  characteristic,  which  was 
experimentally  confirmed  by  Gostick  et  al.  [7],  the  assumption  of 
entirely  hydrophilicity  seems  not  to  be  valid. 

The  mixed  wettability  characteristics  are  caused  by  the  wide 
range  of  wetting  characteristics  of  carbon-based  materials  typically 
used  for  the  PEMFC  GDL,  possible  anomaly  in  the  PTFE  treatment 
and  surface  defects,  impurities  and  aging  of  the  GDL  [21].  Only 
few  researchers  have  considered  the  GDL  as  partially  hydrophobic. 
Weber  et  al.  [23]  proposed  a  composite  contact  angle  as  a  function 
of  the  fraction  of  hydrophilic  pores  in  order  to  take  the  mixed  wet¬ 
tability  of  GDL  into  account,  and  showed  that  there  is  an  optimum 
fraction  value  of  hydrophilic  pores  for  approaching  maximum  lim¬ 
iting  current  and  power.  Nava  et  al.  [12]  determined  the  saturation 


by  averaging  the  values  for  both  the  hydrophobic  and  hydrophilic 
pores  according  to  the  hydrophilic  pore  fraction  and  calculated  the 
capillary  pressure  on  the  dependence  of  the  saturation  value  with 
respect  to  this  fraction.  He  et  al.  [24]  discussed  the  saturations  for 
the  hydrophilicity  and  hydrophobicity  case  respectively  and  also 
averaged  their  values  according  to  the  hydrophilic  pore  fraction. 
Weber  [25]  integrated  analytically  the  separate  hydrophobic  and 
hydrophilic  pore-size  distributions  of  the  diffusion  medium  via  a 
random  cut-and-rejoin  bundle-of-capillaries  model,  and  combined 
their  values  to  get  the  saturation  according  to  the  value  range  of 
contact  angle. 

On  the  view  of  theoretical  research,  the  integral  representation 
method  is  mathematically  more  constitutive  than  others  because 
it  conforms  to  the  pattern  of  liquid  water  movement  through  the 
hydrophilic  and  hydrophobic  pores  rather  than  simply  synthesizes 
according  to  the  hydrophilic  pore  fraction.  For  this  integral  method, 
it  is  crucial  to  determine  the  pore-size  distribution.  Weber  [25] 
assumed  a  series  of  log-normal  distributions  to  fit  the  pore  size  dis¬ 
tribution.  Our  previous  study  [26]  experimentally  provided  insight 
about  the  fractal  characteristic  of  the  GDL  in  PEMFC,  and  the  fractal 
pore  size  distribution  can  be  used  in  the  integral  method. 

The  objective  of  the  present  work  is  to  determine  the  depend¬ 
ence  of  the  capillary  pressure  on  the  saturation  and  the  relative 
permeability  of  gas  and  liquid  water  by  a  fractal  bundle-of- 
hydrophilic-and-hydrophobic-capillaries  model.  And  the  paper  is 
organized  as  follows:  first,  a  fractal  pore-network  model  for  a 
mixed-wetting  carbon  paper  GDL  is  described.  Then,  the  expression 
between  the  capillary  pressure  and  the  saturation  and  the  relative 
permeability  models  of  gas  and  liquid  water  are  presented.  At  last, 
the  parametric  effect  analysis  is  taken. 

2.  Fractal  porous-medium  model  for  the  GDL 

Experimental  investigations  on  permeation  in  porous  media 
have  shown  that  the  channels  through  which  liquids  permeate 
have  fractal  characteristics.  Similarly,  the  gas  transfer  routes  within 
porous  media  may  also  have  fractal  characteristics  and  can  be  rep¬ 
resented  as  random  fractal  curves.  In  PEMFC,  porous  carbon  paper, 
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made  from  graphite  fibers  and  used  for  the  preparation  of  diffusion 
medium,  GDL,  shows  a  fractal  characteristic.  The  gas  and  water  per¬ 
meation  within  the  GDL  can  be  assumed  to  be  similar  to  that  within 
the  tortuous  fractal  parallel  channels  with  different  hydrophilic  and 
hydrophobic  pore  sizes,  depicted  in  Fig.  1. 

The  length  L(  A )  of  the  capillary  pathway  is  related  to  the  capillary 
size  A  (i.e.,  the  pore  diameter)  by  the  following  fractal  relationship 
[27] 

where  L0  is  the  representative  or  linear  length  of  these  capillary 
pathways  towards  the  flowing  direction,  and  Dt  is  the  tortuousity 
fractal  dimension.  Then  L(A)  is 

L(X)  =  L0D<X'~Dt  (2) 


by 


And  the  tortuosity  of  the  capillary  pathway  r  can  be  described 


'=m2=m 


2Dt-2 


(3) 


Another  characteristic  of  porous  media  is  that  the  cumulative 
pore  population  N  in  a  unit  cross  section  may  be  mathematically 
expressed  as  follows  [27]: 

N(d>X)  =  (^p)°P  (4) 

where  Dp  is  the  pore  area  fractal  dimension,  and  A,  Amax  are  the  pore 
size  and  the  maximum  pore  size  of  porous  media,  respectively.  The 
first  derivative  of  Eq.  (10)  with  respect  to  A  is 


-dN  =  DpXfyxX-^P+VdX 


(5) 


where  pL  and  pG  represent  the  liquid  water  and  gas  pressure, 
respectively.  Usually  the  gas  pressure  can  be  considered  as  a  con¬ 
stant  in  the  fuel  cell  [29],  so  the  capillary  pressure  becomes  a  direct 
relation  of  the  liquid  water  pressure. 

For  hydrophilic  case,  the  value  of  Ac  is  the  largest  water- 
filled  pore  diameter,  and  the  contact  angle  is  O°<0<9O°.  As  for 
hydrophobic  case,  Ac  is  the  smallest  water-filled  pore  diameter, 
and  90°  <6><  180°. 

Existence  of  both  hydrophilic  and  hydrophobic  capillaries 
makes  it  difficult  to  investigate  the  wettability  of  GDL.  The  liquid 
saturation  of  diffusion  medium  is  used  for  describing  the  mixed 
wettability,  and  is  defined  as  the  fraction  of  void  space  filled  by 
liquid  water  phase 


spores 

where  vWater  is  the  volume  of  capillaries  or  pores  filled  with  water, 
and  Spores  is  the  total  volume  of  the  tortuous  capillaries  or  pores. 

If  pc  <  0,  i.e.,  at  a  negative  capillary  pressure,  only  the  hydrophilic 
capillaries  with  diameter  less  than  the  critical  diameter  A*j  are 
occupied  by  liquid  water,  then  the  liquid  saturation  s  is 


s  =  - 

=/w 


■/ 


(7T/4)X2L(X)dNhi 


J 


(jt/4  )X2L(X)dN 
Dp  -  X3 


3— Dt— Dp  .3-Dt-Dp 
hi )  min 


y  3-Dt-Dp 
Amax 


y  3-Dt-Dp 
min 


(10) 


As  mentioned  above,  the  porous  diffusion  medium,  GDL,  can  be 
considered  as  a  composite  structure  with  hydrophilic  (usually  car¬ 
bon)  and  hydrophobic  (usually  Teflon)  capillaries,  and  the  number 
of  former  with  diameter  greater  than  A  is  determined  by  [28] 

Nhi(d>^)=/hi(^p)DP  (6) 

where /hi  is  the  fraction  of  hydrophilic  pores,  which  is  determined 
by  the  PTFE  loading. 

The  proposed  fractal  bundle-of-hydrophilic-and-hydrophobic- 
capillaries  model  will  be  used  to  characterize  a  real  carbon  paper. 
Although  this  model  is  not  correct  on  the  local,  pore-scale  level, 
it  is  believed  that  it  remains  valid  at  the  macroscopic,  layer-scale 
level  for  its  enough  description  on  the  real  carbon  paper.  All  the 
following  derivations  are  based  on  this  model. 

3.  Capillary  pressure-liquid  phase  saturation  relationships 
for  diffusion  medium 

Under  the  dry  condition,  the  hydrophilicity  or  hydrophobicity  of 
GDL  has  no  effect  on  the  gas  transport,  but  under  the  wet  condition, 
the  gas  transport  within  GDL  is  highly  related  to  water  transport, 
which  is  in  turn  influenced  by  the  mixed  wettability  of  GDL. 

The  wettability  of  the  GDL  is  related  to  capillary-driven  liquid 
water  behavior,  described  by  the  Young-Laplace  equation 

Xc  =  -^  COS0  (7) 

Pc 

where  Ac  is  the  pore  diameter,  y  is  the  surface  tension  of  water,  0 
is  the  composite  contact  angle  of  water  on  the  solid  surface  of  the 
pore  wall,  and  pc  is  the  capillary  pressure,  calculated  by  [23] 


Pc=Pl-Pg 


(3) 


where  Amin  is  the  minimum  pore  size  of  porous  media. 

If  pc  =  0,  i.e.,  at  a  zero  capillary  pressure,  all  the  hydrophilic  cap¬ 
illaries  are  filled  and  all  the  hydrophobic  capillaries  are  empty,  and 
the  liquid  saturation  s  is 


-  f  lmax(7t/4)X2L(X)dNhi 

£  _  J  Amin _ 

-  f,W(7T/4 )X2L(X)dN 


(11) 


If  pc  >  0,  i.e.,  at  a  positive  capillary  pressure,  all  the  hydrophilic 
capillaries  and  the  hydrophobic  capillaries  with  diameter  greater 
than  the  critical  diameter  A£o  are  occupied,  and  the  liquid  satura¬ 
tion  s  is 


/ 


-r 


s  = 


(7T/4)A2L(A)dNhi  -  /  (7r/4)A2L(A)d(N  -  Nhi) 

A* 

ho 


7 


(7r/4)X2L(X)dN 


^maxt_Dp  -/hiCn''^  -  O  ~  fhiXKo^^ 


(12) 


y  3— Dt— Dp  y  3-Df-Dp 

Amax  “  Amin  ^  ^ 

m*  ^3-Dt-Dp  _  ^  3-Dt-Dp 
h0  min 


y  3— Df— Dp  y 

Amax  -  Amin 


3-Dt-Dp 


From  Eqs.  (1 0)— (1 2),  it  is  easy  to  know  that  the  capillaries  in  the 
GDL  are  occupied  by  more  and  more  liquid  water  with  the  increase 
of  the  capillary  pressure  or  the  liquid  water  pressure,  according  to 
the  following  order  as:  first  the  small  hydrophilic  capillaries,  then 
the  large  hydrophilic  capillaries  followed  by  the  large  hydrophobic 
capillaries,  and  last  the  small  hydrophobic  capillaries. 
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According  to  Eqs.  (10)-(12),  the  liquid  saturation  of  the  GDLcan 
be  described  with  a  piecewise  function 


s  =  < 


/hi 

/hi 


M*  y3-Dt-Dp  _  ,3-Dt-Dp 
lAhU _  min 

y  3-Dt-Dp  ^3-Dt-Dp 


i  -(i  -/w) 


H*  \3-Dt-Dp  _  y  3-Dt-Dp 
^AhJ  Amin 


ho' _ 

\  3—  Dt—  Dp 


-A; 


3-Dt-Dp 


Pc  <  0 

Pc  =  0 

Pc  >  o 


According  to  the  value  range  of  capillary  pressure,  the  total  flow 
rate  of  water  Qw  is 


q(A)dNhi 


(13)  Qw  —  < 


■f 

7 

Mmax  M 

-  /  q(X)dNhi  -  / 

<  ^min 


''-min 
/>Amax 

q(A)dNh: 

/^/nax 

q(l)dNhi-  /  q(X)d(N  -  Nhi)  pc  >  0 

^min  "  “ho 


Pc  <0 

Pc  =  0  (19) 


And  then  the  relation  between  the  capillary  pressure  and  the 
liquid  saturation  can  be  obtained  by  substituting  Eq.  (7)  into  Eq. 

(13) 


s  =  < 


(-((4K  cos  6>hi)/Pc))3  Dt  Dp  -  A3-Dt-Dp 


n  3-Dt-Dp  n  3-Dt-Dp 


i  -d  -Ai) 


(-((4y  cos  0ho)/pc)) 


3— Df— Dp  ^3— D^—  Dp 


^max'"Dp  -  A3TDt"Dp 
max  min 


Pc  <  0 
Pc  =  0 

Pc  >  0 


(14) 


where  0hi  is  the  contact  angle  for  the  hydrophilic  phase  and  0ho  is 
the  contact  angle  for  the  hydrophobic  coating. 

From  Eqs.  (13)  and  (14),  it  can  be  seen  that  the  liquid  saturation 
of  the  diffusion  medium  is  a  continuous  function  of  the  capillary 
pressure,  and  ifpc<0,  since  A.* j  <  Amax,  then  s</hi;  if pc  >0,  because 
A,*o  <  ^max,  s  >fhi,  and  Eq.  (14)  can  also  be  rewritten  as 


Using  Darcy’s  law,  we  can  obtain  the  permeability  for  saturated 
porous  diffusion  medium 


/xL0Q 
A  pA 


(20) 


and  the  permeability  for  liquid  water 


K\v  — 


Mw^oOw 

ApwA 


(21) 


where  /xw  is  the  viscosity  of  water,  A pw  is  the  pressure  gradient  of 
water,  and  A  is  the  representative  area. 

Substituting  Eqs.  (2),  (5),  (17)  and  (18)  into  Eq.  (20),  we  have 


rl—  Dt 


I<  = 


n 

U'nAma’s, 


128  A  3-Dp+Dt 


(  y3-Dp+Dt  _  y3-Dp+Dt\ 
vAmax  Amin  > 


(22) 


And  according  to  Eqs.  (2),  (5),  (6),  (17),  (19)  and  (21),  we  obtain 


s  =  < 


/hi 


(-(4k  cos  0hi/Pc))3"D,"Dp  -A3"“'"Dp 


'-XT 


s  </hi 


(-(4y  cos  0ho/Pc))3~Dt-Dp  -  ^~Dp 
1  (1  /hi)  _  o_n._n„  .  S  >/hi 


K 


3— Df—  Dp  ^3-Dt-Dp 


(15) 


Because  there  is  no  liquid  transport  at  the  residual  liquid  satura¬ 
tion,  s r,  an  effective  liquid  saturation,  se,  which  can  be  interpreted 
as  the  fraction  of  transportable  liquid,  has  to  be  taken  into  account 
as  follows: 


I<w  —  i 


n 

r  1— Dt 

/hi£+max 

128 

A 

3-Dp+D, 

71 

T  1— Dt 

/hi£+max 

128 

A 

3  —  Dp  +  Dt 

71 

rl  -Dt 

Dp  max 

128  A  3-Dp+Dt 

-(1  -/hi)(A.*0)3-6p+Dt] 


[(KJ3~DP+Dt  -Am7nP+Dtl 

(A.^“°P+Dt  -  X3~DP+Dt ) 

v  max  min  J 

[&+Dt  -/hi^3m7nP+Dt 


Pc  <  0 
Pc  =  0 


Pc  >  0 


(23) 


Combining  Eqs.  (22),  (23)  and  the  definition  of  the  water  relative 
permeability  yields 


se  — 


S  -  Sr 
1  —  Sr 


(16) 


4.  Fractal  model  on  relative  permeability  in  diffusion 
medium 

4.1.  Fractal  model  on  water  relative  permeability  in  diffusion 
medium 


The  water  relative  permeability  is  of  great  importance  for 
species  simulation  in  the  GDL.  The  flow  rate  q(X)  through  a  sin¬ 
gle  tortuous  capillary  with  a  diameter,  A,  is  given  by  modifying  the 
well  known  Hagen-Poiseulle  equation  [30] 


q(X)  = 


71  A p  A4 

128  L(X)  7/ 


(17) 


I<w 

~K 


,  3-Dp+Dt 


Pc  <  0 
Pc  =  0 
Pc  >  0 


(24) 


From  Eq.  (24),  it  can  be  seen  that  the  water  relative  perme¬ 
ability  is  a  continuous  function  of  the  capillary  pressure,  and  if 
pc  <  0,  since  A.JL  <  A,max,  then  /<rw</hil  if  Pc>0,  because  A.min  < 
A,*o  <  ^max,/hi  <  k rw  <1,  and  Eq.  (24)  can  be  reduced  as 


/hi 


/ y  *  ^3-Dp+Dt  y  3  Dp+Dt 

(A.J  -Amin 


/<rw  —  < 


/hi  + 


y  3— Dp+Dt  _  y 
Amax  7  Amin 

(i  -/w) 


3— Dp  +Dt 


^3maxP+Dt-(^0)3“Dp+Dt 


>  3— Dp+Dt  y  3— Dp+Dt 

A™;„ 


krw  <  /hi 


krw  >  /hi 


(25) 


where  A  P  is  the  pressure  gradient,  and  /x  is  the  viscosity  of  the  fluid. 
Then  the  total  volumetric  flow  rate  Qcan  be  obtained  by  integrating 
the  individual  flow  rate  q(A)  over  the  entire  range  of  pore  sizes 
from  the  minimum  pore  size  A,min  to  the  maximum  pore  size  Amax 
in  a  unit  cell  according  to  the  fractal  bundle-of-hydrophilic-and- 
hydrophobic-capillaries  model,  as  follows: 

/A  max 

Q  =  -/  q(X)dN  (18) 

^  ^-min 


According  to  Eq.  (25),  the  water  relative  permeability  is  also  a 
continuous  function  of  the  hydrophilic  pore  fraction /hi. 

4.2.  Fractal  model  on  the  relative  permeability  of  gas  in  diffusion 
medium 

In  the  simulation  of  PEMFC,  the  gas  permeability  diffusion  in 
GDL  also  needs  to  be  taken  into  account. 

At  a  negative  capillary  pressure,  only  the  hydrophilic  capillaries 
with  diameter  less  than  the  critical  diameter  X*  are  occupied  by 
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Table  1 

Values  for  model  parameters  and  physical  properties  of  GDL  material  tested. 


Parameter 

Value 

Hydrophilic  pore  fraction 

/hi 

0.35 

Contact  angle  for  water  on  graphite  [23] 

Ohi 

80° 

Contact  angle  for  water  on  PTFE  [23] 

Oho 

110° 

Residual  liquid  saturation  [25] 

ST 

0.05 

Water  surface  tension  [23] 

Y 

0.12398-0.00017393TN  m-1 

Maximum  pore  diameter  [24] 

^-max 

1.0661  x  10“4m 

Minimum  pore  diameter 

A,  min 

2.5  x  10-6  m 

Pore  area  dimension  [31  ] 

DP 

1.9669 

Tortuosity  dimension  [31  ] 

Dt 

1.1447 

Operation  temperature 

T 

333  K 

Absolute  permeability  [31] 

I< 

8  x  10-12  m2 

Porosity  [31] 

s 

0.78 

liquid  water,  gas  can  pass  through  the  rest  capillaries,  and  then  the 
total  flow  rate  of  gas  Qg  is 


Qg  -  - 


rXmax  rx  w 

/  q{X)dN  +  /  q{X)dNh[ 

0  Amin  "  ^min 


(26) 


At  a  zero  capillary  pressure,  all  the  hydrophilic  capillaries  are 
filled  with  water,  gas  can  pass  through  all  the  hydrophobic  capil¬ 
laries,  and  the  total  flow  rate  of  gas  Qg  is 

/*A  max 

Qg  =  -/  q(A)d(N-Nhi)  (27) 


For  a  positive  capillary  pressure,  all  the  hydrophilic  capillaries 
and  the  hydrophobic  capillaries  with  diameter  greater  than  the  crit¬ 
ical  diameter  A*o  are  occupied  by  water  and  the  total  flow  rate  of 
gas  Qg  is 


Qg  =  -  / 


q(X)d(N-Nhi) 


(28) 


Combining  Eqs.  (2),  (5),  (6),  (26),  (27)  and  (28),  the  gas  perme¬ 
ability  can  be  gained  by 


Kg=  < 


n 

r\-Dt 

Lo 

Dpi 

128 

A 

3 -L 

-fM(K+ 

71 

.1-1+ 

L0 

(1-. 

128 

A 

3- 

71 

.1-Dt 

L0 

(1-. 

128 

A 

3- 

-  [A,3-Dp+Dt 

L^-max 

'-(i  -/hiU^>p+1><] 
i  Dp 


Pc  <  0 
Pc  =  0 


(29) 


And  the  relative  permeability  of  gas  is 


k  -K*-l 


1  “/hi  + 
1  -/hi 
(1  “/hi) 


/hi[^maxP+Dt-(^i)3"Dp+D‘] 


h  3— Dp+Dt  ^  3— Dp+Dt 
Amax  -  Amin 

f  y  *  ^3— Dp+Dt  Tl  3— Dp+Dt 

(AhJ  Amin 


ho^ _ 

•\  3 -Dp+Dt  y  3— Dp+Dt 
Amax  _  Amin 


reduced  as 

fhi[^maxP+Dt  -(^hi)3-Pp+Pt] 


Pc  <  0 
Pc  =  0  (30) 
Pc  >  0 


,  r  y  3— Dp+Dt 

Amax  —  Amin 
n3— Dp  +Dt 
o' 


y  3— Dp+Dt 
Amax 


ltv _ 

3— Dp  +Dt 


^  j— ivp-rivt 

,  35)ip+Dt 

Amin 


krg  >  1  -/hi 


5.  Results  and  discussion 

From  the  above-mentioned  derivations,  we  can  see  that  the  pc-s 
correlation  is  significant  to  the  proposed  models.  As  far  as  the  pc-s 


Table  2 

Fitting  parameters  for  various  empirical  relations. 


Parameter 

Value 

Arbi 

0.8040  kPa 

Ab2 

0.7655  kPa 

n 

15.2839 

m 

0.0678 

V 

1.0016 

correlation  is  concerned,  a  mention  should  be  made  concerning 
hysteresis,  which  is  common  in  the  porous  media  literature.  Sev¬ 
eral  researchers  [25,31,32]  have  observed  the  hysteretic  wetting 
between  liquid  and  gas  intrusion  into  the  GDL,  which  demonstrates 
preferential  wetting  of  the  pores,  the  importance  of  history,  and 
that  the  GDL  pores  exhibit  intermediate  wettability  [25  ].  It  is  largely 
due  to  “ink-bottle”  effects  and  hysteresis  in  the  contact  angle  [31  ]. 
Here  we  have  derived  the  pc-s  correlation  according  to  the  process 
of  liquid  intrusion,  without  regard  to  that  of  gas  intrusion. 

From  the  above  fractal  models  of  saturation  (see  Eqs.  (14)  and 
(15))  and  relative  permeability  of  water  and  gas  (see  Eqs.  (24),  (25), 
(30)  and  (31 )),  it  can  be  seen  that  both  the  liquid  water  saturation 
and  the  relative  permeability  are  determined  not  only  by  capil¬ 
lary  pressure,  contact  angles  and  water  surface  tension  but  also  by 
the  GDL  structure  parameters  (/hi,  ^max,  A,min,  Dp  and  Dt).  For  the 
following  theoretical  analysis,  a  set  of  parameters  and  operating 
conditions  has  been  specified  as  the  base  case,  which  comes  from 
the  literatures  [23-25,33],  as  listed  in  Table  1. 

5.1.  Predicted  saturation 


The  saturation  depends  on  the  capillary  pressure,  which  is  a 
quantitative  measure  of  the  liquid  water  transport  in  the  GDL  and 
is  related  to  the  diffusive  transport  of  liquid  water  since  the  liquid 
water  diffusion  coefficient  is  defined  in  terms  of  the  capillary  pres¬ 
sure  [20].  Several  models  are  commonly  used  to  fitpc-s  curves.  The 
most  widely  used  one  is  the  Leverett  J-function  model  [20] 

Pc  =  X  cos  6»(|:)  /(s)  (32) 


where  s  is  the  air-filled  porosity,  I< is  the  absolute  permeability,  and 
the  term/(s),  called  the  Leverett  J-function,  represents  the  dimen¬ 
sionless  capillary  pressure  as  a  function  of  liquid  saturation  and  is 
given  as 

fn=  f  1.417(1  -s)-  2.120(1  -s)2  +  1.263(1  -s)3  if  0°  <  0  <  90°  /oon 

/U  \  1.417s -2.120s2  +  1.263s3  if  90°<6><180°  ^  ] 

Other  common  models  are  the  van  Genuchten  model  (Eq.  (34)) 
and  the  Brooks-Corey  model  (Eq.  (35))  [7]. 


Pc  >  0 


(34) 
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s=(f^)  '■  Pc>Pcb2  (35) 

where  pc bi ,  pCb2»  n,  m,  v  are  fitting  parameters. 

The  present  fractal  model  is  compared  with  the  above  empirical 
models  and  the  experimental  data  complied  from  the  literature 
[31],  as  shown  in  Fig.  2.  These  experimental  data  are  obtained 
from  the  process  of  liquid  intrusion,  which  corresponds  with 
the  assumption  of  our  models.  The  fitting  parameters  of  the  van 
Genuchten  model  or  the  Brooks-Corey  model  can  be  determined 


from  a  best  fit  of  the  model  to  experimental  data,  as  given  in  Table  2. 
The  best  fit  is  determined  by  the  lowest  root  sum  of  squared  error 
between  the  experimental  data  and  the  capillary  model.  And  the 
value  of  hydrophilic  pore  fraction  in  the  fractal  model  is  assumed  to 
0.35  for  the  hydrophobic  tendency  of  the  experimental  sample.  As 
can  be  seen,  the  results  of  our  fractal  model  are  in  better  agreement 
with  the  experimental  data  than  other  empirical  models.  Hereinto, 
the  Leverett  J-function  model  overestimates  the  saturation,  and  the 
curve  slope  of  this  function  is  only  close  to  that  of  experimental 


Fig.  3.  Liquid  water  saturations  as  function  of  capillary  pressure  for  different  hydrophilic  pore  fraction. 
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Fig.  4.  Effect  of  tortuosity  fractal  dimension  Dt  on  liquid  water  saturation. 


data  during  a  small  range.  As  for  the  van  Genuchten  model  and  the 
Brooks-Corey  model,  data  points  below  zero  capillary  pressure  and 
Pcb2  are  excluded  respectively,  since  these  models  cannot  describe 
the  behavior  in  these  regions.  Although  these  empirical  models 
agree  well  with  part  of  experimental  data  for  the  hydrophilic  case, 
they  show  a  bad  consistent  with  the  experimental  data  in  the  case 
of  mixed  wettability,  depicted  in  Fig.  2.  Furthermore,  a  drawback 
of  all  these  empirical  models  is  that  they  do  not  reveal  any  direct 
relationship  between  pc-s  curve  and  the  GDL  structure. 


From  Fig.  2,  it  also  can  be  seen  that  the  saturation  varies 
slowly  with  negative  capillary  pressure  while  varies  greatly  with 
positive  capillary  pressure  because  the  sample  with  a  small 
hydrophilic  pore  fraction  includes  few  hydrophilic-like  pores  corre¬ 
sponding  to  negative  capillary  pressure  and  many  hydrophobic-like 
pores  corresponding  to  positive  capillary.  Moreover,  the  satura¬ 
tion  increases  as  capillary  pressure  increases  from  negative  value 
to  positive  one. 

In  order  to  investigate  the  effect  of  PTFE  on  the  pc-s  curve,  differ¬ 
ent  values  of  hydrophilic  pore  fraction  were  chosen,  as  illustrated  in 


Fig.  5.  Effect  of  pore  area  fractal  dimension  Dp  on  liquid  water  saturation. 


Y.  Shi  et  al.  /  Journal  of  Power  Sources  209  (2012)  130-140 


137 


Fig.  3.  As  can  be  seen,  the  value  of  saturation  increases  greatly  as  the 
hydrophilic  pore  fraction/hi  increases  from  0.2  to  0.8.  This  increase 
should  correspond  to  more  hydrophilic  pores  and  larger  porosity 
on  lower  PTFE  loadings.  At  the  same  time,  the  value  of  capillary 
pressure  increases  as  the  hydrophilic  pore  fraction  /hi  decreases, 
which  means  that  capillary  pressure  increases  with  the  addition  of 
PTFE.  This  increase  is  caused  by  the  decrease  of  the  total  pore  space 
available  for  the  water  to  occupy  and  the  increase  in  the  amount  of 
hydrophobic  pores. 


In  addition,  from  Fig.  3,  it  also  can  be  seen  that  at  a  hydrophilic 
pore  fraction  lower  than  0.5,  i.e.,  for  a  nearly  hydrophobic  case,  the 
capillary  pressure  increases  slowly  with  increasing  s  at  low  liquid 
saturations  but  quickly  when  saturation  is  higher.  It  is  due  to  that 
few  hydrophilic  capillaries  exist  and  liquid  water  prefers  to  fill  the 
large  hydrophobic  capillaries.  While  at  a  hydrophilic  pore  fraction 
higher  than  0.5,  i.e.,  for  a  nearly  hydrophilic  case,  the  opposite  result 
is  obtained  because  a  lot  of  hydrophilic  capillaries  exist  and  the 
small  hydrophilic  capillaries  are  first  to  be  filled. 


Fig.  7.  Water  relative  permeabilities  as  function  of  water  saturation  for  different  hydrophilic  pore  fraction. 
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Fig.  3  also  shows  that  for  the  mixed  wettability  GDL  the  variation 
ofpc-s  curve  exhibits  a  trend  of  a  slow  increase  at  lower  and  higher 
capillary  pressure  but  a  fast  raise  around  zero  capillary  pressure. 
The  occupation  of  small  hydrophilic  capillaries  at  the  beginning  of 
condensation  leads  to  the  slow  increase  of  pc-s  curve  at  negative 
capillary  pressure.  After  this,  the  fast  raise  of  pc-s  curve  around 
zero  capillary  pressure  is  attributed  to  the  liquid  water  filling  in 
the  large  hydrophilic  capillaries  as  well  as  the  large  hydrophobic 
capillaries.  And  then  the  slow  increase  of  pc-s  curve  at  positive 


capillary  pressure  is  contributed  by  the  occupation  of  the  rest  small 
hydrophobic  capillaries. 

The  effects  of  tortuosity  and  pore  area  fractal  dimensions  onpc-s 
curve  were  investigated,  as  shown  in  Figs.  4  and  5,  respectively.  It 
can  be  seen  that  for  a  hydrophilic  case  the  value  of  capillary  pressure 
decreases  with  the  increase  of  fractal  dimensions  at  a  given  satura¬ 
tion  and  the  value  of  saturation  increases  as  the  fractal  dimensions 
increase  at  a  given  capillary  pressure  while  for  a  hydrophobic  case 
the  opposite  result  is  observed.  From  Eq.  (1)  and  Eq.  (4),  we  can 
know  that  a  larger  tortuosity  and  pore  area  fractal  dimensions 


Fig.  9.  Effect  of  tortuosity  fractal  dimension  Dt  on  water  relative  permeability. 
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Fig.  10.  Effect  of  pore  area  fractal  dimension  Dp  on  water  relative  permeability. 


correspond  to  a  longer  capillary  length  and  an  increasing  capil¬ 
lary  number,  respectively.  The  above  results  indicate  that  higher 
fractal  dimensions  result  in  greater  increments  of  the  number  and 
length  for  smaller  capillaries  and  a  lower  fraction  of  pore  space 
contributed  by  the  large  capillaries.  Therefore,  for  the  hydrophilic 
GDL  with  higher  fractal  dimensions,  more  small  capillaries  will  be 
occupied  by  liquid  water  at  a  given  capillary  pressure,  leading  to 
a  higher  saturation.  While  for  the  hydrophobic  case,  the  opposite 
conclusion  can  be  drawn,  this  is  caused  by  the  first  occupation  of 
large  capillaries. 

5.2.  Predicted  relative  permeability 

The  fractal  prediction  on  gas  and  water  relative  permeability 
for /hi  =  0.2  and  0.8  is  presented  in  Fig.  6.  It  is  seen  that  the  relative 
permeability  of  the  liquid  phase  is  zero  until  liquid  water  satura¬ 
tion  exceeds  the  threshold  value  at  residual  liquid  saturation.  Fig.  6 
also  shows  that  for  the  GDL  of  the  mixed  wettability,  the  behavior  of 
water  relative  permeability  versus  saturation  is  controlled  by  three 
distinct  regimes,  i.e.,  a  slow  increase  regime,  a  fast  raise  regime  and 
again  a  slow  increase  regime,  whereas  that  of  gas  relative  perme¬ 
ability  versus  saturation  is  controlled  by  three  other  regimes,  i.e.,  a 
slow  decrease  regimes,  a  fast  drop  regime  and  again  a  slow  reduc¬ 
tion  regime.  The  explanation  for  this  result  is  similar  to  that  for  the 
variation  ofpc-s  curve  in  the  mixed-wettability  case. 

Figs.  7  and  8  show  that  the  variations  in  relative  permeabil¬ 
ity  with  the  saturation  for  different  hydrophilic  pore  fractions.  The 
water  relative  permeability  increases  and  the  gas  relative  perme¬ 
ability  decreases  with  the  reduction  of  hydrophilic  pore  fraction. 
It  is  also  seen  that  water  relative  permeability  of  the  hydropho¬ 
bic  case  is  much  larger  than  that  of  the  hydrophilic  case;  while  for 
gas  relative  permeability  the  opposite  result  is  obtained.  Therefore, 
for  the  goal  to  remove  liquid  water,  a  hydrophobic  GDL  is  more 
suitable  than  the  hydrophilic  one  for  its  high  water  relative  perme¬ 
ability.  While  for  the  goal  to  maintain  a  certain  quantity  of  water, 
a  hydrophilic  GDL  is  better  for  its  low  water  relative  permeability. 

Figs.  9  and  10  show  that  the  water  relative  permeability  varies 
with  the  saturation  for  different  values  of  tortuosity  and  pore  area 


fractal  dimensions,  respectively.  The  results  indicate  that  the  water 
relative  permeability  of  the  hydrophilic  case  decreases  and  that 
of  the  hydrophobic  case  increases  as  the  two  fractal  dimensions 
increase,  which  has  been  observed  by  the  study  of  Fie  et  al.  [24]. 

6.  Conclusions 

The  purpose  of  this  study  is  to  propose  theoretical  models  to 
determine  the  saturation,  water  and  gas  relative  permeability  of 
the  GDL  in  PEMFCs.  To  this  end,  the  relationships  between  the  sat¬ 
uration,  relative  permeability  and  the  microstructure  of  the  GDL 
are  considered  using  fractal  methods,  respectively.  To  validate  the 
proposed  fractal  model  of  saturation  versus  capillary  pressure,  we 
compare  it  to  experimental  data  and  three  empirical  models  that 
are  widely  used  in  the  research  of  pc-s  curve.  It  can  be  concluded 
that  in  the  case  of  mixed  wettability,  our  fractal  model  is  better  for  a 
good  agreement  with  the  experimental  data,  whereas  these  empir¬ 
ical  models  are  not  suitable  for  the  GDL  of  mixed  wettability  due 
to  obvious  deviation  between  their  results  and  the  experimental 
data.  After  this,  the  effects  of  PTFE,  tortuosity  and  pore  area  frac¬ 
tal  dimensions  on  pc-s  curve  are  investigated.  It  is  seen  that  the 
saturation  is  positively  correlated  with  the  two  fractal  dimensions 
for  a  hydrophilic  case  and  the  hydrophilic  pore  fraction,  whereas 
it  is  negatively  correlated  with  the  two  fractal  dimensions  for  a 
hydrophobic  case.  And  the  capillary  pressure  is  negatively  corre¬ 
lated  with  the  two  fractal  dimensions  for  a  hydrophilic  GDL  and  the 
hydrophilic  pore  fraction,  whereas  it  is  positively  corrected  with 
the  two  fractal  dimensions  for  a  hydrophobic  GDL.  These  results 
indicate  that  higher  fractal  dimensions  corresponds  to  greater 
increments  of  the  number  and  length  for  smaller  capillaries  and 
a  larger  fraction  of  pore  space  contributed  by  the  small  capillaries. 
Furthermore,  the  prediction  on  gas  and  water  relative  permeability 
and  the  analysis  on  parametric  effect  are  taken  through  the  pro¬ 
posed  models  of  water  and  gas  relative  permeability.  Water  relative 
permeability  increases  with  the  increases  in  the  two  fractal  dimen¬ 
sions  for  the  hydrophobic  GDL  and  liquid  water  saturation,  whereas 
it  decreases  with  the  increases  in  these  fractal  dimensions  for  the 
hydrophilic  GDL  and  hydrophilic  pore  fraction  at  a  given  saturation. 
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As  for  gas  phase,  gas  relative  permeability  is  positively  correlated 
with  hydrophilic  pore  fraction,  while  it  is  negatively  correlated  with 
liquid  water  saturation. 

The  proposed  models  provide  a  bridge  between  microscale  and 
macroscale  models  of  PEMFCs.  In  future  studies,  we  will  leverage 
our  previous  work  on  the  properties  of  the  GDL  to  perform  multi¬ 
scale  and  two-phase  modeling  of  PEMFCs  and  thus  facilitate  the 
design  of  the  GDL  and  of  PEMFCs. 
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